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1. Introduction 



After the release of the BaBar data [[[]] for the pion-photon transition form factor, F^o™., our 
community is in deep shock, as the conventional approach to the gold-plated exclusive process 
at very high Euclidean momenta Q, based on 1. factorization and 2. (leading-twist) pQCD evolu- 
tion |§-|T2|], seems to be invalid [13, 14]. Indeed, the quantity Q 2 Fjpyy*(Q 2 ) goes visibly above 
the famous Brodsky-Lepage limit, 



y ' N c Jo 



(1.1) 



at momenta Q 2 > 15 GeV 2 . 

Several ideas, abandoning assumptions 1. and 2., have been proposed to solve the "BaBar 
problem". Radyushkin [|T^] and Polyakov [jlfj advocated that the possible nonvanishing of the 
PDA at the end points (as found by the present authors in chiral quark models at the low-energy 



quark-model scale Q17Q), together with essentially switched-off evolution and regulated quark prop- 
agators, is capable of reproducing the data in the CLEO and BaBar domain. In this approach 
Fjfiyyt ~ log(2 2 / i u 2 )/2 2 , with the log indicating the breaking of factorization (note that the same 
asymptotics follows in the Spectral Quark Model (SQM), cf. Eq. (14.1) of [Q). Dorokhov [|^- 
21] proposed the use of the fixed-mass chiral quark model to evaluate the triangle diagram of 
Fig. [I], which is capable of reproducing the data at high Q 2 , however the fitted value of the con- 
stituent quark mass is very low, M ~ 135 MeV. In this the form Fjfi yy , ~ [log(<2 2 //i 2 )] 2 /<2 2 . A 
possible need of higher-twist terms has been brought up in p2[]. The calculation of [p3|] in the 



nonlocal chiral quark model inspired by the instanton-liquid model of QCD produced the result in 
agreement with the data at lower values of Q 2 and complying to the limit (JO]). On the other hand, 

Z5D in the nonlocal model and found 



Dorokhov [ ]24p considered a modified pion-meson vertex 
agreement in the whole available data range, albeit again with a very low constituent quark mass. 
Asymptotically, in this approach fjjo™* ~ log(2 2 //i 2 )/2 2 . The influence of nonperturbative glu- 
onic components of the pion has been considered in Concerns on the validity of factorization 
were discussed in [[2^]. It was also shown that the pronounced growth of Q 2 F 7l 0yy> between 10 and 



40 GeV cannot be explained via higher-order pQCD corrections at the NNLO level [12]. Further 



developments may be found in [28 - 31 ] 
This talk is based on [f32|l. 




Figure 1: The quark-model diagrams used to evaluate the pion transition form factor (left) and the Zq — > n y 
decay (right). The crossed diagrams not shown. 
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2. General constraints 



Consider the process of the left part of Fig. [I] with the general kinematics 



q\ = ~ l -^Q 2 , <& = - l —^Q 2 ' -1 <A< 1. (2.1) 



The chiral anomaly [33, 34] fixes F^o^fiQ = 0,A) = -^rj-- 

On the other, the high-<2 2 behavior of F n 0yy> is formally limited by the Terazawa-West (TW) [35 



37Q unitarity bounds, recently brought up by Dorokhov pip and further elaborated in [|32||. The 
derivation of the TW bounds uses the Schwarz inequality involving sums of the matrix elements 
(0|/u(0)|n) and (7i a (q) \Jn(0)\n), entering the vacuum polarization and the partem distribution func- 
tion (PDF) of the pion. Schematically, one has |(?r|//|0)| < (|(0|77|0)||(7r|77|7r)|) 1/2 . Then one 
derives the first bound, 

ImF 7r o rr (^ 2 ) = ^(l/7?) (TWI) (2.2) 



holding at time-like momenta, q 2 > 4m 2 . If there are no polynomial terms in the real part of F n ayy* 
(which is an assumption), then |^V> rr *(g 2 )| = ^(l/^/q 2 ). A dispersion relation yields [37] the 
second bound, 

\F n o rr (Q 2 ) \ = 6{\jQ) (TWH) (2.3) 
valid for all momenta, also large space-like momenta Q. The constant in the bound may be explic- 



itly given Q37| ] in terms of the photon spectral density and the pion structure function, 



\F^ (fi ) | < — ^— J q dx^j ^ (2.4) 

where H(s) = s/(l6n 3 'aQ ED )a e + e -^i ia( i Tons (s) andll(oo) = \/(l2n 2 )J^ i ej, withe,- denting the quark 
charges. With the SMRS [ |38| ] and GRV [ [3^ ] parameterizations for F\ we obtain (for Q 2 in the range 
10-40 GeV 2 ) 

l^o rr (e 2 )|<^ (LO), <^p (NLO), (2.5) 

where LO and NLO refer to Leading Order and Next-to-Leading Order in the QCD evolution, 
respectively. Thus the TW II bound is "inefficient" in the present experimental range, as it goes an 
order of magnitude above the BaBar data. 

Finally, another interesting experimental bound comes from the rare Z — > 7i°y decay pCp, 
which probes the time-like value q 2 = M z . The process is shown in the right part of Fig. [T[ Since 
only the vector coupling of the Z° boson to the quark contributes, 

(2.6) 



The experimental limit given by the Particle Data Group igj, T(Z -> 7T°y) < 5 x 10" 5 r tot (Z°) = 
10.25 x 10~ 5 GeV, implies 

\F z o^o r (M 2 )/F z o^o r (0)\<0.n. (2.7) 
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3. Subtracted dispersion relation and violation of TW II 

The assumption of the absence of the polynomial terms, necessary for TW II to hold, is equiv- 
alent to the validity of the unsubtracted dispersion relation for F^oyy. Clearly, pQCD with factor- 
ization leads to F„0y*y vanishing at Q — > °°, however, since now these assumptions are questioned 
(see the Introduction), one may consider the situation where subtraction constants are necessary in 
the dispersion relation for iyy r []32|]. Note that even if the form factor vanishes at infinity, one can 
write a subtracted relation 

1 f 4Al t lmF(s) 1 r t ImF(s) 

F(t)-F(0) = - —ds + - —ds. (3.1) 

% Js s s-t K Jak 2 s s-t 

If A is large (A 2 > Q 2 ), the second term is very slowly varying with Q 2 and mimics a constant. In 
particular, for ImF(s) ~ l/y^it behaves as 1/A + &(l/Q). Thus the appearance of the constant 
term need not be taken as a fundamental problem, as it may represent the unknown high-energy 
data. We will analyze below a quantitative lower bound for a possible high energy mass scale. 

The bottom line is that one may well consider the case where F n a y *y{Q 2 ) does not vanish at 
Q 2 — > oo. Below we will provide an explicit field-theoretic example where this is the case (the 
Georgi-Manohar (GM) model p2|]). Motivated by this we will also incorporate a constant term in 
the fits of the BaBar data. 



4. Mini-review of chiral quark models 

Before discussing the GM model let us briefly review the chiral quark models and, in partic- 
ular, their predictions for the soft matrix elements entering the (factorized) high-energy processes. 
The purpose is to convince the reader that these models, based on chiral symmetry breaking as 
the key dynamical ingredient, are very successful for numerous processes involving the pions and 
photons. Chiral quark models are cast in a covariant Lagrangian form and in general exhibit no 
factorization. They carry relatively few parameters, traded for f n , m n , . . . The large-A^. limit is 
implicit, however, confinement is absent and one needs to be careful not to open the qq production 
threshold. This limits the applicability range, however, no problems arise with space-like external 
momenta, as in the present analysis of the BaBar puzzle. Finally, predictions of chiral quark mod- 



els hold at a low-energy quark-model scale [43] and QCD evolution is necessary to reach higher 
scales of experiments or lattice simulations. 



We start with the valence PDF of the pion, where the NJL model gives [44] 



q(x) = 1 (4.1) 



(note that the same constant PDF for two constituents follows from AdS/CFT approach of [45]). 



Arguments based on the momentum sum rule [43] determine the renormalization scale where 



Eq. (gTj) holds: /i ~ 320 MeV. Then at the scale \i=2 GeV the valence quarks carry 47% of 
the momentum of the pion, as requested by experimentally motivated parameterizations. The value 
of the coupling constant at the quark-model scale is a(jlio) jit = 0.68. After the evolution the NJL 
prediction agrees remarkably well with the experimental data (left panel of Fig. |2|). 
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xq(x) 




Figure 2: Left: The NJL prediction for the valence PDF of the pion evolved to the experimental scale 
ji = 4 GeV (band). The data points come from the analysis of the Drell-Yan data from the E615 experiment 



6h. The dashed line shows the reanalysis of the data from [47]. Right: The NJL prediction for the PDA 



of the pion evolved to the lattice scale n = 0.5GeV (band) compared to the transverse lattice data [ 48 ]. The 
width of the bands indicates the uncertainty in the quark-model scale jUq. 



In the NJL model, the PDA of the pion is also constant at the quark-model scale [|l 

0to = 1 



(4-2) 



(note that this result is different from the AdS/CFT prediction of ty(x) ~ \/x(l — x), it is also 
far from the asymptotic form 6x(l —x)). The PDA ( |4.2[ ) does not vanish at the end points, which is 
the focal point of Jl5| , 16]. However, the LO ERBL evolution makes 0(a) vanish at the end-points, 
with the form (j)(x) ~ x^r/h i°g[«(«>)/a(jO] near x = [g [43]]. Thus maintaining the nonvanishing 
at the end points requires switching off the QCD evolution fll5[]. 

Finally, we mention the NJL results for the gravitational form factor of the pion, discussed 



in detail in [49]. The quark-model relation holds for the radii related to the spin-2 gravitational 



form factor and the charge form factor, {r 2 )s = \{r 2 )v- Therefore matter is more concentrated 



than charge, in agreement with the recent AdS/CFT results Q50|]. Note that there is no contradiction 
between having a constant vertex function and a finite pion size. This issue is tightly linked with 



the fulfillment of the electromagnetic and chiral Ward identities, as discussed, e.g., in [51]. 



5. Violating TW II 

We now recall a model which reproduces the results of the previous section, but violates TW II: 



the Georgi-Manohar (GM) model [12]. It is obtained from the NJL model by carrying out the chiral 
rotation (which is innocuous) and then introducing gA of the quark which may be different from 



unity [52]. The Lagrangian of the GM model is 



L = q (id+gQAYs-M^Jq+t-Tr (d^d^U) + WZW, 



where WZW stands for the Wess-Zumino-Witten term, while 



— (u'duU — udnU^) 



,*?/(2/) 



u 



(5.1) 



(5.2) 
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Figure 3: Left: predictions of the one-state IVMD model (line). Right: predictions of various Regge 
models: the Veneziano-Dominguez model b = 1.81 (dashed line), model with the first pole separated and 
fixed b = 1.5 (dotted line), and the subtracted Regge model (solid line). The dots, squares, and diamonds 
correspond to the CELLO pB]], CLEO and BaBar |[IJ data, correspondingly. 



The resulting pion-photon transition form factor is 
1 



■ + 



\% 2 f n 4% 2 f K 



g Q A T / 7s 1 / 7s 2Ml f l dx 

' M [G(Q 2 )-\], G(Q) 



Q 1 



f l dx 
■ / —log 

Jo x 



2 1 



(5.3) 



As we can see, the anomaly is satisfied, but for gf / 1 no vanishing at Q 2 — > oo occurs, 

- {Q 2y \-8 Q A , gfMMlQg(QVM 2 )] 2 | _ 



In the Spectral Quark Model JT^ ] one has instead 



Q 2 



(5.4) 



G(e 2 ) = ~ 



2m 2 p 



m~ mi, + <2 



m 2 +Q 2 ' 2 



> ^7C° 



1 



?m 2 log; 



rt 



■ + 



A% 2 f % \2% 2 f n Q 2 



(5.5) 



With gQ = 1 this model fulfills the result of [15] with the mass scale m p . 

Unfortunately, precise chiral-quark-model fits of F K o rr based on above formulas in the whole 
Q 2 range are not satisfactory, so in the following we proceed to more general analyses based on 
Incomplete Vector Meson Dominance (IVMD) [53, 54]. 



6. Incomplete VMD and Regge models 

As shown in the previous Section, it is possible to formulate a field-theoretic model con- 
sistent with all formal requirements which violates TW II, leading to an asymptotically constant 
Fjfiyyt (<2 2 )- With this finding in mind we explore the idea of IVMD and the Regge models. Regge 
models with (infinitely many) tree-level meson and glueball exchanges are a realization of the 
large-A^ limit which, unlike the current quark models, incorporate confinement and comply with 
quark-hadron duality by generating meromorphic two-point functions for color singlet currents. 
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The simplest IVMD model which comes to mind has just one state saturating the subuacted 



dispersion relation for F n 0yy,.. Then 



71° ff 



(Q 2 



Q 



'M 2 + Q 2 



(6.1) 



The fit to the combined CELLO CLEO [Q, and BaBar [g] data yields c = 0.986(2), M v = 
748(14) MeV, with x 2 /DOF = 0.7, thus (1 -c) is significantly different from although nu- 
merically small. On the other hand, the fit to the CLEO data only gives c = 0.998(18), Mv = 
777(44) MeV, £ 2 /DOF = 0.54, with (1 - c) compatible with 0, or complete VMD. The result of 
the fit to the combined data is shown in the left panel of Fig. |3[ displaying a remarkable agreement 
in the whole experimentally available range. The radius squared 



F n o yr (<2) 



F^ rf {Q)d& 71 yr 



Q 2 =0 



(6.2) 



becomes b n 



A/r 



1.76(7) GeV 2 , compared to the PDG value b n = (1.76±0.22)GeV~ 2 . The 



constrain from the rare Z° decay is satisfied comfortably, as \F z ^ n 0y(M z ) / F^^oJO)] = 0.014(2), 



an order of magnitude less than 0. 17 of Eq. ( |2.7[ ). 

Next, we consider radial Regge models, with M 2 = M v + an, recalling that the large-/V c QCD 
involves tree-level diagrams with infinitely many states, including the radial excitations. A particu- 
lar realization, the Veneziano-Dominguez model [57, ^], allows to control the asymptotic behavior 
of the form factor with a single parameter, b: 

1 



Frf) yr (t) 



4tt 2 A 
1 



fb{t) 



b + n) 



B(b-l,^) 



y T(2 

hr(n + i)r(2 -b)M% 



i 



(6.3) 



n=0 



where fb(0) = 1, fb(t = —Q 2 ) ~ (Q 2 ) l b - The same approach was used in [ |59| ] to successfully 
describe the pion charge form factor, also at large Q 2 . Note that b < 1.5 complies to TW II, 
while b > 1.5 violates the bound. In the right panel of Fig. ^ we show the results of several 
Regge models: the Veneziano-Dominguez model with the fitted value b = l.Sl (dashed line), the 
Veneziano-Dominguez model with the first pole separated and fixed b = 1.5 (dotted line), and the 



subtracted Regge model (solid line). The details can be found in Q32p . We note that the models 
differ only at very large values of Q 2 , where the present experimental uncertainties are large. 

Hence, it is possible to explain the BaBar data within the Regge approach, both in the TW II 
violating and TW II conserving versions. 



Along the lines of the discussion around Eq. (p.l|), another adventurous scenario can be pro- 
posed, which lies between the one mass state case and the infinitely many uniformly spaced squared 
mass states. We can then interpret the constant c for the IVMD ansatz of Eq. ( |6.1[ ) as corresponding 
to a high-mass state, M#, such that for M 2 , <C Q 2 <C Mjj one has a Q 2 -independent behavior. This 
can be represented by the simple functional form 



l 



4tt 2 A 



My 



'M 2 + Q 2 



M% 



M 2 +Q 2 



(6.4) 
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which resembles the IVMD ansatz of Eq. dO) in the range M v <C Q 2 <C M#. A direct fit to the 
joint CLEO and BaBar data, with fixed M v = 748 MeV, yields c = 0.085(2) and the lowest bound 
(within the 95% confidence level) for the high-mass state is M# ~ 10 GeV. 



7. Photon momentum asymmetry 

The last issue we wish to discuss in this talk is the influence of the photon momentum asym- 
metry parameter, A. In the BaBar kinematic setup — q\ < 0.6 GeV 2 and —q\ > 3 GeV 2 , therefore 



\A\ 



2 2 



qj + ql 



~ 0.9 -0.97 / 1. 



(7.1) 



This turn out to be significant for precision fits and the optimum values of their physical parameters. 
In particular, the single-state IMVD model becomes 



7T' yy* 



(-Q 2 



l 



4tt 2 A 



1-c 1 



4M V 



4M v +4M v Q 2 + (l-A 2 )Q 4 



(7.2) 



The fits to the combined data yield for A = 1, 0.975, and 0.95, respectively, c = 0.986, 0.978, 0.974, 
My = 748, 754, 768 MeV. In particular, the value of My increases toward m p as A decreases. 



8. Conclusions 



These are our main conclusions: 



• The second Terazawa-West bound, following from assumptions for the real part of the pion 
transition form factor, need not be fulfilled in field-theoretic approaches. An explicit counter- 
example with the Georgi-Manohar model with g& / 1 , where the chiral anomaly is fulfilled 
but F K 0yy, tends to a constant at t — > — oo. If this feature holds in QCD, it opens a possibil- 
ity of solving the BaBar puzzle within models incorporating the incomplete vector-meson 
dominance. 

• The coefficient in the TW II is large, extending the bound an order of magnitude above the 
data. Thus, even if it holds in the absence of polynomial contributions to the real part, it is 
completely ineffective for the presently-available momentum range. 

• An additional constraint on the models in the time-like region of momenta follows from the 
rare Z — > 7t°j decay. We use this bound in our considerations. For the considered models it 
is comfortably satisfied. 

• The simplest model realizing the incomplete vector-meson dominance, including a single 
vector-meson state, is capable of reproducing the data for F n 0yy,- in the whole available ex- 
perimental range, < Q 2 < 35 GeV 2 . 

• Within the Regge approach with infinitely many radially excited states, the data can be fitted 
both satisfying or violating the TW II bound. 
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• The precise numerical fits are sensitive to the photon momentum asymmetry parameter, A, 
which affect the optimum values of the physical parameters, such as the vector meson mass. 
Since A is not strictly 1 in the experimental setup, the effects of kinematic cuts should be 
considered in precision analyses. 
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